Quantum Fluctuations in the Chirped Pendulum 
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Anharmonic oscillators, such as the pendulum, 
are widely used for precision measurement [lj 
and to model nonlinear phenomena [2j. In any 
physical detector, fluctuations — such as thermal 
or quantum mechanical noise — can excite random 
motion in the oscillator, ultimately imposing a 
bound on measurement sensitivity. In systems 
where equilibrium is established with the noisy 
environment, for example most electrical or op- 
tical amplifiers, noise-induced broadening scales 
with the intensity of fluctuations, which for ther- 
mal noise is proportional to temperature (T). But 
how does noise affect an out of equilibrium os- 
cillator where the motion is varied faster than 
energy is exchanged with the environment? We 
create such a scenario by applying fast, frequency 
chirped voltage pulses to a nonlinear supercon- 
ducting resonator where the ring down time is 
longer than the pulse duration. Under these con- 
ditions, the circuit oscillates with either small or 
large amplitude depending on whether the drive 
voltage is below or above a critical value— a phe- 
nomenon known as autoresonance |3j. This sharp 
threshold is significant in planetary dynamics |4j 
and plasmas [5j, enables the excitation of parti- 
cles in cyclotron accelerators |6j , and may even be 
used to detect the state of a quantum two level 
system [7J. For the first time, we observe a noise- 
induced broadening of the autoresonant thresh- 
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old that is ex T e ^ where T e ff = T in the classical 
regime and hw/2kB in the quantum regime. The 
results imply that fluctuations only determine the 
initial conditions of such a non-equilibrium oscil- 
lator and do not effect its time evolution. 

Our oscillator is a superconducting electrical circuit 
based on a Josephson tunnel junction formed by two thin 
aluminum films separated by an oxide barrier. In re- 
sponse to a periodic current drive /(t), the junction volt- 
age V(t) oscillates according to the Josephson relations, 
which describe the motion of a particle with coordinate S 
in a sinusoidal potential U(5) = [hIo/(2e)] cos(^), where 
Io is the junction critical current, 6 is the phase difference 
across the junction, and e is the electron charge. This sys- 
tem is analogous to a mechanical pendulum with angular 
coordinate — 5 (Fig. la and b). The response of such an 
anharmonic oscillator is illustrated in Figure lc; as the 
drive strength is increased the resonant response shifts 
to lower frequency and ultimately exhibits bistability. In 



this regime, if the oscillator is excited with a constant 
frequency drive of increasing amplitude, as depicted by 
the light grey arrow in Figure lc, it abruptly switches 
from a low to a high amplitude oscillation state. This 
effect has been used to realize a threshold amplifier [8] to 
readout the state of quantum bits [9HTT] . In the presence 
of fluctuations, the switching from one oscillation state to 
the other is modeled as an equilibrium Arrhenius process 
involving escape over an effective energy barrier [I2j [13] . 
Consequently, the sensitivity of this amplifier improves 
with decreasing temperature until the quantum regime 
is reached [13] . 

The bistable regime may also be accessed by driving 
the anharmonic oscillator with a fixed amplitude, chirped 
drive as depicted by the black arrows in Figure lc. For a 
given chirp rate a and oscillator quality factor Q, if the 
drive amplitude is above a critical voltage Vc a 3 / 4 , the 
oscillator phase locks with the drive signal and climbs to 
the high amplitude oscillation state [3 . Although known 
in plasma and accelerator physics since the 1950s, au- 
toresonance has only recently been observed in an electri- 
cal circuit, where it was noted that this process could be 
used for sensitive amplification [7]. Unlike the amplitude 
driven case, autoresonant excitation is by construction 
a non-equilibrium process, and in this letter we analyze 
the finite width of the threshold and the role of quantum 
fluctuations in determining the ultimate sensitivity of a 
measurement device which utilizes this effect. 

The oscillator circuit consists of a tunnel junction em- 
bedded in the center of a linear cavity (see Fig. Id), re- 
alized using a capacitively isolated section of coplanar 
waveguide transmission line that forms a 6 GHz Fabry- 
Perot resonator. This architecture allows us to engineer 
the electrical impedance shunting the junction and hence 
tune the oscillator quality factor (Q) and provides a sim- 
ple means to couple microwave frequency signals into and 
out of the oscillator. The equation of motion for the em- 
bedded system at large excitation differs from that of a 
physical pendulum; however, the autoresonant dynamics 
are unchanged. We excite the circuit with a frequency 
chirp, created using a voltage controlled oscillator modu- 
lated by an arbitrary waveform generator. This signal is 
sent via heavily filtered coaxial lines to the chip, which is 
cooled to temperatures as low as 15 mK in a dilution re- 
frigerator. The transmitted voltage is amplified at 2.3 K 
and demodulated using quadrature mixers at room tem- 
perature to obtain the magnitude relative to the drive. 

We first demonstrate threshold behavior when the os- 
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FIG. 1: The pendulum. The resonant response of an anharmonic oscillator, as realized by a pendulum (a) or a capacitively 
shunted Josephson tunnel junction (b), becomes increasingly asymmetric (c) as the drive strength is increased. Above a critical 
excitation the oscillator exhibits two metastable states, as indicated by black circles, with low and high amplitude of oscillation. 
These states may be excited by ramping the drive amplitude at a fixed frequency (grey dashed arrow) or by chirping the drive 
at fixed power (black dashed arrow), d, Schematic of the oscillator circuit consisting of a superconducting resonator embedded 
with a Josephson junction. The frequency modulated drive signal is generated using a voltage controlled oscillator (VCO) and 
detected by means of a cryogenic amplifier and room temperature electronics (not shown). 



cillator is driven with a phase continuous voltage pulse 
that starts at 6.075 GHz and decreases linearly in time 
to 5.775 GHz. The linear resonance frequency of the os- 
cillator is uj/(2tt) = 5.987 GHz with Q = 8200. Thus, 
both the start and stop frequencies are many linewidths 
away from resonance. The transmitted voltage (V) de- 
veloped in response to the drive is plotted in Figure 2a as 
a function of drive frequency and amplitude. The data 
are acquired by averaging 5,000 frequency sweeps at a 
fixed amplitude and then stepping its value over a 10 
nV interval. As the drive becomes commensurate with 
the resonance frequency, oscillation amplitude begins to 
build. When the drive amplitude is weak, we observe 
that after a slight excitation the oscillator relaxes back 
to very small oscillations. On the other hand, for stronger 
drive, the oscillator builds up energy, behavior that is in- 
dicative of autoresonant phase-locking of the oscillations 
to the drive. 

A constant frequency slice of the data in Figure 2a, 
shown in the left panel, shows the threshold for au- 
toresonance at V c ~ 27 nV. Here, the average junction 
voltage varies between a low level (Vi) corresponding to 
the unlocked events and a high level (Vh) correspond- 
ing to locked events. Within the threshold, in a given 
frequency sweep the oscillator either locks into autoreso- 
nance or does not. Because we average many sweeps, we 
extract the locking probability from this measurement, 
P = (V - Vi)/(V h - Vi). The locking probability ver- 
sus drive power is shown for three different chirp rates in 
Figure 2b. The threshold location, V c , where P = 1/2, 
is plotted versus chirp rate in Figure 2c, and agrees with 
theory [7]. 

We now turn to the width of the autoresonant thresh- 



old. In the absence of fluctuations, the capture probabil- 
ity would be an infinitely sharp step at V c as indicated in 
Figures 2d and 2e. At the start of a chirp sequence, the 
resonator has an initial excitation that is the result of in- 
teractions with the fluctuating bath up to that time. The 
threshold for a given initial amplitude Aq is thus either 
augmented or diminished depending on initial the phase 
mismatch, A</>, between the drive signal and the oscilla- 
tor motion, V c = V c — kAq cos A(/>, where n = 0.245 is a 
prefactor determined by numerical simulations. A finite 
threshold width is obtained by integrating the locking 
probability (see Fig. 2d) over all values of A0 and over a 
thermal distribution of initial amplitudes at temperature 
T eff , where k B T eS = ^-^-t'/tfr] - The threshold 
width, defined as the inverse slope of the locking proba- 
bility at P = 1/2 is then AV C = 2n\/2n La/c^Te^, where 
L is the total inductance of the resonator and ks is Boltz- 
mann's constant [T5] . 

To confirm the validity of this description when ksT <C 
Huj, we compute the dynamics in the presence of quan- 
tum fluctuations. The state of an oscillator can be de- 
scribed by its position and momentum, which for our 
electrical circuit correspond to the dimensionless charge 
(q) and current (j), respectively. In Figure 3, we plot 
the Wigner quasi-probability distribution that results af- 
ter the ground state — a Gaussian centered about zero — is 
evolved using a voltage excitation 20 % greater than V c . 
When the anharmonicity is large compared to dissipa- 
tion, high excitation is achieved by climbing the ladder 
of accessible states, as opposed to a continuous classi- 
cal evolution [16 . Panels a-c correspond to decreasing 
anharmonicity with energy level shifts 100, 49, and 25 
times larger than our experimental parameters; panel d 



3 




FIG. 2: Autoresonance in the chirped pendulum, a, The average transmitted voltage (V) versus drive frequency and 
drive voltage. The drive was chirped at a rate of 50.6 MHz//xs. The dashed line indicates the fixed frequency used for the slice 
(left panel), which shows the threshold for autoresonance to occur at V c — 27 nV. b, The locking probability ascertained near 
the autoresonant threshold for chirp rates of 42.1 (blue), 50.6(black), and 63.2 (red) MHz//is. c, The threshold location, V c 
varies with the chirp rate and agrees with theoretical predictions 7 . d, Cartoon of the locking/non-locking regions for finite 
amplitude (solid line) and zero (dashed line) fluctuation versus the phase mismatch between the drive and the oscillator, e, 
Averaging over the initial amplitudes and phases of fluctuations broadens the transition (solid line) from the zero fluctuation 
(dashed line) prediction. 
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FIG. 3: Calculated Wigner distributions. Calculated 
Wigner distributions (f(q,j,r = 4215)) for the quantum (a- 
c) and classical nonlinear oscillator (d) as a function of the 
dimensionless charge q and current j. a-c correspond to non- 
linear energy level shifts that are respectively 100, 49, and 25 
times larger than in the experiment. The classical result in d 
shows that approximately 90 % of the distribution is in the 
phase-locked state, as indicated by the dashed line. The tail 
outside this region relaxes to low amplitude if dissipation is 
included. Negative values in the classical calculation result 
from the finite accuracy of the simulation. 



is the classical result with the majority of the population 
in the phase locked state, as indicated by the dashed 
line. The data in Figure 3c show nearly the same dis- 
tribution as that of the classical calculation. For higher 
anharmonicity, the phase locked state remains qualita- 
tively unchanged, but the population in the low ampli- 
tude state exhibits a characteristic quantum interference 
pattern indicating the participation of a relatively small 
number of discrete levels. Our sample parameters, how- 
ever, correspond to an oscillator with 25 times weaker 
anharmonicity than Figure 3c, and we thus expect the 
classical model for the transition width to be valid at all 
temperatures. 

By varying the physical temperature of the dilution re- 
frigerator we measured the transition width as a function 
of temperature. To compare to theory we plot the square 
of the transition width, scaled to temperature units as 
AV 2 j (SKirLakB) , versus temperature in Figure 4. For 
temperatures above 200 mK, we observe a clear T 1 / 2 
dependence of the width. At the lowest temperatures, 
saturation is observed corresponding to T e R = 144 mK, 
which is precisely fkj/(2/c#). The solid line indicates our 
theoretical prediction for the width with no adjustable 
parameters. To verify the quantum origin of this satu- 
ration, we plot the threshold width data for a 1.6 GHz 
resonator (Figure 4 inset). Here, we observe a T 1 / 2 scal- 
ing over the entire temperature range, with a suggestive 
flattening of the transition width at the lowest tempera- 
tures, where fc#T c± hw. 

Keeping the sample at T = 15 mK, we also injected 
white noise into the resonator to simulate a thermal bath. 
The results of this experiment are shown as the blue open 
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FIG. 4: The autoresonance threshold width. The 

square of the threshold width, scaled to temperature units 
as AV 2 /(8k,7t Laths), and averaged across chirp rates rang- 
ing from 50 to 500 MHz//xs versus system temperature (red 
circles), saturates to the value predicted for quantum fluc- 
tuations. The data are in agreement with our predictions, 
represented as an effective temperature, T e ff, and shown as a 
solid black line. White noise injected into the system causes 
the threshold width to increase oc T 1 ^ 2 (blue open circles) 
with an offset given by the zero point quantum fluctuations 
in the resonator. Inset: the threshold width for a junction 
embedded resonator at 1.6 GHz exhibits oc T 1 ^ 2 scaling down 
to T = 40 inK. 



circles in Figure 4. The observed threshold width scales 
as T 1 / 2 and has an offset corresponding to the half quan- 
tum of zero point motion of the oscillator, corroborat- 
ing again a quantum noise broadened threshold at the 
lowest temperature. Thus, quantum mechanics in our 
system enters only through the initial fluctuations, char- 
acterized by T e ff , and the system otherwise can be treated 
classically [T7] . 

Our experiment is an example of the interaction of an 
oscillator out of equilibrium with a noise bath. The sim- 
ple temperature scaling of the threshold width can poten- 
tially be harnessed for robust noise thermometry. Fur- 
thermore, the quantum saturation of this width, which 
we note can be tuned by adjusting a, ultimately sets the 
resolution of a digital detector based on autoresonance. 
Such a detector can be used for the readout of quantum 
bits. In contrast to abrupt switching between metastable 
states [18], autoresonance involves continuous evolution 
to either the locked or unlocked state, raising interesting 
questions about the measurement backaction. 

Methods The oscillator was fabricated from Al de- 
posited on a high-resistivity Si wafer. At low excitation 



power, a total quality factor of Q = 8230 resulted from 
internal losses characterized by Qi nt = 17200 and cou- 
pling to the 50 Q environment given by Q ex t — 15800. 
The critical power for bifurcation was measured to be 
P c = —123 dBm, corresponding to a junction critical 
current Io = 1.8 fiA and a total inductance L = 2.3 nH. 
At drive powers near P c we observed a decrease in the 
total quality factor of the oscillator to Q = 6480. 

The chirped pulse was created by driving a VCO oper- 
ating near 8 GHz with an arbitrary waveform generator 
and mixing it with a tone at 2.312 GHz. The transmitted 
power from the sample was mixed with the output of the 
VCO to form a fixed frequency signal at 2.312 GHz. This 
signal was demodulated with a local oscillator at 2.412 
GHz, and the I and Q quadratures were digitized at 500 
MS/s. 

We modeled our oscillator via a weakly nonlinear 
Hamiltonian, H = j 2 /2 + q 2 /2 — (3q 4 /4 + eq cos <\>&, on 
the phase space of dimensionless charge, g, normalized to 
Qo = jo/w an d current, j, normalized to jo — ^Jk^T^jL. 
The anharmonicity is given by f3 = ($>quj 2 Qq)/(6LIq), 
where $0 = h/2e is the flux quantum. The dimension- 
less drive voltage e = Vd/ '(Lq^uj 2 ) and the phase of the 
drive <pd = ojt — at 2 / 2. We numerically solved the di- 
mensionless quantum Liouville equation for the Wigner 
function, f(q,j,r), where r = cut: 



[(<? - Pq 3 ) +e cos <f>d 



df 7 2 Az a 3 / 



dj 



4 df 



(1) 



The parameter 7 = two / '(fce^eff) allows us to solve both 
the classical 7—^0 and quantum 7—^2 limits of the sys- 
tem. We solved the Liouville equation using a standard 
pseudospectral method [19] for the parameters in our ex- 



periment, but for fj = n x 3.55 x 10 



0.0246/n, 



and a = 10~ 6 u 2 Hz 2 , where n = {10,7,5} correspond 
to nonlinear energy level shifts that are 100, 49 and 25 
times that of the experiment, and drive strengths that 
are 20% above the autoresonant threshold. The classi- 
cal evolution was approximated by solving the Liouville 
equation with 7 = 10 -4 . 
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